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Abstract 

We consider contributions of R-parity conserving softly broken su- 
persymmetry (SUSY) to neutrinoless double beta (Oi^/3/3) decay via the 
(B-L)-violating sneutrino mass term. The latter is a generic ingredient 
of any weak-scale SUSY model with a Major ana neutrino mass. The 
new R-parity conserving SUSY contributions to Qv(3(3 are realized at 
the level of box diagrams. We derive the effective Lagrangian describ- 
ing the SUSY-box mechanism of Oi^/^/J-decay and the corresponding 
nuclear matrix elements. The 1-loop sneutrino contribution to the Ma- 
jorana neutrino mass is also derived. 

Given the data on the Oi^/3/3-decay half-life of ^^Ge and the neu- 
trino mass we obtain constraints on the (B-L)-violating sneutrino mass. 
These constraints leave room for accelerator searches for certain mani- 
festations of the 2nd and 3rd generation (B-L)-violating sneutrino mass 
term, but are most probably too tight for first generation (B-L)-violating 
sneutrino masses to be searched for directly. 
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1 Introduction 



Neutrinoless double beta (Oz//3/3) decay 0| is a unique example of a nuclear 
process which allows one to probe for lepton number violation. Given the fact 
that the standard model conserves L at the classical level, 0^(3(3 decay can 
proceed only via non-SM interactions. Therefore - taking into account the 
stringent experimental limits - 0^(3/3 decay is extremely sensitive to physics 
beyond the standard model. 

The simplest source of lepton number violation directly leading to 0^(3(3 
decay is a finite Majorana neutrino rest mass. It violates lepton number by two 
units AL = 2 precisely what is necessary for Oi^/3/3-decay. The corresponding 
contribution is described by the tree-level diagram of the 4th order in the 
weak coupling constant with one Majorana neutrino propagator as shown in 
Fig.l(a). 

Supersymmetric extensions of the SM bring in new sources of lepton num- 
ber violation and, as a result, new mechanisms of 0z//3/3-decay. Even with the 
minimal matter and Higgs field content there are renormalizable L and B vi- 
olating terms in the superpotential and in the sector of soft-supersymmetry 
breaking interactions which are not forbidden by gauge symmetry. These terms 
violate not only L and B but also R-parity defined as Rp = (^— i)3'S+l+25 ^j^]-^ 
5* being a particle spin. Models containing such interaction terms are usually 
referred to as SUSY models with explicit _Rp-breaking {Rp can also be bro- 
ken spontaneously via a non-zero vacuum expectation value of the sneutrino 
field < z> >^ |].) 

Formerly it was widely believed that supersymmetry can contribute to 
Oi>f3(3 decay only if Rp is broken The corresponding mechanisms have 



been comprehensively studied in the last few years 0|-[|TT]. In particular, it 
was shown that the current experimental limits from non-observation of Qv(3(3- 
decay sets upper bounds on certain If,p Yukawa coupling constants which are 
more stringent |^, ^ than previously known from various accelerator and non- 
accelerator experiments. Moreover, they turned out to be more stringent than 
those from some forthcoming accelerator experiments. This conclusion has put 
Oi//5/5-decay forward as an interesting probe of supersymmetry. 

Although there are no compelling theoretical arguments for R-parity con- 
servation, there exist a number of well-known phenomenological drawbacks for 
supersymmetric models in which i?p is violated. Maybe the most serious one 
is the instability of the lightest SUSY particle. As a result, the supersym- 
metric solution for the dark matter problem is lost unless Yukawa coupling 
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constants A become unnaturally small, typically A < 10~^^. 

In view of this and other problems for ^pSUSY there arises a natural ques- 
tion whether i?p-violation is an inevitable condition for SUSY to contribute 
to Oz//3/5-decay. The present paper addresses this question. We will demon- 
strate that there is a non-trivial R-parity conserving SUSY contribution to 
Oi//3/5-decay. In Fig. 1(b) we present an example of a diagram associated with 
the lowest order contribution to Oz//3/3-decay within the R-parity conserving 
minimal supersymmetric standard model (MSSM) with a Majorana neutrino 
mass m'^M- This particular example gives an explicit answer to the above ques- 
tion: R-parity violation is not a necessary condition for a SUSY contribution 
to Oz//5/5-decay. 

It is also obvious that this SUSY contribution is strongly suppressed com- 
pared to the non-SUSY diagram in Fig. 1(a). This is because it is of higher 
order in perturbation theory, contains heavy sparticles in intermediate states 
and receives a typical suppression due to the loop integration. Moreover, this 
diagram is proportional to the very small factor mlj/pp where pp ~ 80MeV 
is the nucleon Fermi momentum. The latter is also true for the simplest non- 
SUSY diagram in Fig. 1(a). The reason is common for both diagrams. In fact, 
both the SM and the MSSM interactions conserve lepton number L. Therefore, 
the only source for AL = 2 violation, necessary for Oz//5/3-decay to proceed, 
is the Majorana neutrino mass term. If m'^.j = lepton number would be a 
conserved quantity and the Oz//3/3-decay amplitude should vanish, = 0. 

Such a behavior corresponds to Rouisp ~ 'm^^/pp in the limit of small m\j. 

Thus, the diagram in Fig. 1(b), given its very small contribution to Qv(3(3- 
decay, provides just a principal demonstration of the fact that Oi^/5/3-decay 
can be triggered by R-parity conserving supersymmetry. No practical conse- 
quences can be obtained from this new diagram in the sense of establishing 
new constraints either on SUSY parameters or on from non-observation 
of Oz//3/3-decay. 

However, we will show that there are other R-parity conserving SUSY con- 
tributions via the lepton number violating sneutrino mass term. As was shown 
r2[ the Majorana neutrino mass, the (B-L)-violating sneutrino mass and 
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the Oz//3/3-decay amplitude are generically connected to each other. Namely, 
non-vanishing of one of these three quantities implies non-zero values of the re- 
maining two. Thus, the Ot//9/3-amplitude should always contain a contribution 
corresponding to the (B-L)-violating sneutrino mass term. 

We will study this contribution and extract constraints on the (B-L)- 
violating sneutrino mass term from the current experimental data on 0z//9/3- 
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decay of Ge. 

This paper is organized as follows. In section 2 we give a short account 
on the structure of neutrino-sneutrino mass terms and the theorem which es- 
tablishes the above mentioned relations between the neutrino and sneutrino 
masses and Oz//3/3-decay. Sect. 3 is devoted to some general properties of pos- 
sible SUSY contributions to 0z//9/9-decay. In this section we specify the box 
diagrams describing the i?p- conserving SUSY contribution. Our approach to 
the derivation of the corresponding Oi//5/?-transition operators and nuclear ma- 
trix elements is outlined in sect. 4. Sect. 5 deals with Oi//3/3-decay constraints 
on the (B-L)-violating sneutrino mass. In section 6 we calculate the sneu- 
trino contribution to the Majorana neutrino mass and derive then limits from 
experimental data on neutrino masses. We then close with a short summary. 



Structure of the Neutrino-Sneutrino Mass 
Terms 



As shown in [|12| , the self-consistent form of the neutrino and sneutrino mass 
terms is 

Qlass = -\i.TAi^^L + h.C.) - ^{mljULyi + h.C.) - ml^luL- (1) 

where z/ = i/^ is a Majorana field. The first two terms violate the global 
(B-L) symmetry while the last one respects it. The first term is a Majorana 
mass term of the neutrino. We call the second term a "Majorana" -like mass, 
while the third one is referred to as a "Dirac"-like sneutrino mass term. This 
reflects an analogy with Majorana and Dirac mass terms for neutrinos. In the 
presence of the right handed neutrino field the Dirac neutrino mass term 
TTT'oi^L^R + ^r^l) could also be included in Eq. (|l]) but it is not required 
by the self-consistency arguments. Note that is not a positively defined 
parameter. 

Eq. (|l]) is a generic consequence of weak-scale softly broken supersymmetry 
and does not depend on the specific mechanism of mass generation in the low- 
energy theory. For the sake of simplicity and without any loss of generality we 
ignore possible neutrino mixing. 



The low-energy theorem proven in ref. |T2[ relates the following three (B- 
L)-violating quantities: the neutrino Majorana mass m'^j, the "Majorana"-like 
sneutrino mass ttim and the amplitude of Oz//5/5-decay Rqu^p- Here we shortly 
describe the proof of this theorem. 



4 



It is relatively easy to see that if at least one of the quantities is non- 
zero the two others are generated in higher orders of perturbation theory as 
demonstrated in Fig. 2, where only dominant diagrams are shown. Inter- 
nal lines in these diagrams are neutralinos Xiy gluinos g, charginos x^y selec- 
tron e, u-squark u and sneutrino u. The latter is to be identified with the 
(B-L)-violating "Majorana" propagator proportional to m\.j. The sneutrino 



"Majorana" propagator was explicitly derived in ref. [|T2[ and will be given 
below. 

The various diagrams lead to relations among the three (B-L)-violating 
observables, which we write down schematically 

Zi ^ ^ fljj ■ Zj -\- (2) 

Here, Zi can stand for Zi = mlj, rh\.j, Rou/3/3- The coefficients ttij correspond to 
contributions of the diagrams in Fig.2(a)-(f) so that i,j = a, b, c, d, e, f. Terms 
Ai represent any other possible contributions. The explicit form of Qij and Ai 
is not essential in the following. Important is only the presence of a correlation 
between m^^, m\^, Roupis, expressed by Eq. (H). 

Now we are going to prove that if = then Zi^ = = (the same 
will be true for any permutation). On the basis of Fig. 2 and Eqs. (|]) one 
can expect such properties of the set of observables Zj. Indeed Zi^ = in the 
left-hand side of Eq. strongly disfavors zj^ ^ and z^^ ^ 0, because it 
requires either all the three terms in the right-hand sides to vanish or their 
net cancelation. The latter is "unnatural". Even if such a cancelation would 
be done by hand, using (unnatural) fine-tuning of certain parameters, in some 
specific order of perturbation theory, it would be spoiled again in higher orders 
of perturbation theory. The cancelation of all terms in the right-hand side of 
Eqs. (0) in all orders of perturbation theory could only be guaranteed by a 
special unbroken symmetry. Let us envisage this possibility in details. 

The effective Lagrangian of a generic model of weak scale softly broken 
supersymmetry contains after electro-weak symmetry breaking the following 
terms ||TB| 



^ = -V2gei- uiXii^L - ■ CLXi i^L ~ get ■ i^Lxt(^L+ (3) 
+ -^{t^lYcl + ULYdL)W; + g ■ x^Y{0^,PL + O^^Pr)xIW- 

+ ... + h.c. 

Dots denote other terms which are not essential for our further considerations. 
Here, ul and represent scalar superpartners of the left-handed neutrino ul 
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and electron fields. The chargino xf ^ind neutralino Xi superpositions 
of the gaugino and the higgsino fields. The contents of these superpositions 
depend on the model. Note that the neutralino is a Majorana field Xi = Xi- 
The explicit form of the coefficients ej, ef and Ofj^ is also unessential. For the 
case of the MSSM one can find them, for instance in |13|. Eq. (j^) is a general 



consequence of the underlying weak scale softly broken supersymmetry and 
the spontaneously broken electro-weak gauge symmetry. 

The Lagrangian (^) does not posses any continuous symmetry having non- 
trivial (B-L) transformation properties. Recall, that U{1)b-l is assumed to 
be broken since we admit (B-L)-violating mass terms in Eq. (|1]). However, 
there might be an appropriate unbroken discrete symmetry. Let us specify this 
discrete symmetry group by the following field transformations 

u -> r]^u, z/ r]ci), cl Veei, ei Veei, (4) 
W+ Vw^^^ Xi Vx^X^, X^ Vx+X^, Ql VqQL- 

Here rji are phase factors. Since the Lagrangian (j^) is assumed to be invariant 
under these transformations one obtains the following relations 

vtVuVx. = 1' VeVx+V*u = 1. ••• (5) 
VeVwVt = 1. VwVx+V*x. = 1' VdVwVu = 1 

Dots denote other relations which are not essential here. The complete set 
of these equations defines the admissible discrete symmetry group of the La- 
grangian in Eq. (|^). 

Let us find the transformation property of the operator structure respon- 
sible for Oz//3/5-decay under this group. At the quark level Oi//3/3-decay implies 
the transition dd — >• uuee, described by the effective operator 

Co./3/3 = a, ■ uT^pd ■ urf )rf ■ erf ^e^ (6) 

where ctj are numerical constants, rf ■* are certain combinations of Dirac 7 
matrices. The Oz//3/5-decay amplitude -R01//3/3 is related to the matrix element 
of this operator 

i?o./3/3~<2e-(A,Z + 2)|Oo./3/3|(A^) > (7) 

where {A, Z) is a nucleus with the atomic weight A and the total charge Z. 
The operator in Eq. @ transforms under the group @ as follows 

OquPP rilupp^Qupp (8) 
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with 



Voupp = VdVuVl (9) 
Solving Eqs. @, (||), one finds 

vl = vl = ivlppY- (10) 

This relation proves the statements 1,2. To see this we note that the observable 
quantity Zi = (m^j, Th\j, Roupis) is forbidden by this symmetry if the corre- 
sponding discrete group factor is non-trivial, i.e. rjf ^ \. Contrary, if rjf = 1, 
this quantity is not protected by the symmetry and appears in higher orders of 
perturbation theory, even if it is not included at the tree-level. Relation (|TUp 
claims that if one of the Zi is forbidden then the two others are also forbidden 
and, vice versa, if one of them is not forbidden they are all not forbidden. 

This completes the proof of the theorem relating the neutrino Majorana 
mass m'^,j, the "Majorana" -like sneutrino mass rfiM and the amplitude of Qv(3(3- 
decay Rou/3i3- The proven theorem can be considered as a supersymmetric 
generalization of the well-known theorem relating only neutrino Majorana 
mass and the neutrinoless double beta decay amplitude. 

Let us show that the "Dirac"-like (B-L) conserving sneutrino mass "rfiD 
should also be present in the theory to ensure the stability of the vacuum 
state. Towards this end consider the last two terms of Eq. (|l]) which we 
denote as C^^g^ and use the real field representation for the complex scalar 
sneutrino field 



i>= (z>i+iz>2)/V2, (11) 
where i>i 2 are real fields. Then 

ri/ l/~2~~ ,7, \ ~2~*~ 1~2~2 1~2~2 ^io\ 



where 



"^1,2 = "^D 



± I^mI (13) 



Assume the vacuum state is stable. Then 2 — i-^- ""^d — I^mI' otherwise 
the vacuum is unstable and subsequent spontaneous symmetry breaking occurs 
via non-zero vacuum expectation values of the sneutrino fields < Vi >^ 0. The 
broken symmetry in this case is the R-parity. It is a discrete symmetry defined 
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as Rp = (— i)3-B+i+25^ where S, B and L are the spin, the baryon and the 
lepton quantum number. 

Therefore, as indicated at the beginning of this section the self-consistent 
structure of the mass terms of the neutrino-sneutrino sector is given by Eq. (||) . 
The mass parameter friM gives a measure of sneutrino-antisneutrino mixing 
{v — u*). This (B-L)-violating effect is an evident manifestation of the 2nd 
term in Eq. (|l|). On the other hand, a finite rfiM gives rise to sphtting the 
complex scalar field z/ = (z/i + iv-i)! ^2 into two real mass eigenfields v\^2 with 
the masses = "^d I^mI- According to the above definition, v\ is the 
CP-even state while is the CP-odd one. 

Let us write down the explicit form of the above mentioned (B-L)-violating 
"Majorana" propagator A^^ for the sneutrino []T2[ which is necessary for our 



subsequent considerations. It can be derived by the use of the real field rep- 
resentation as in Eq. (|T^. For comparison we also give the (B-L)-conserving 
"Dirac" Af* sneutrino propagator, 

A^(x -y) = i< 0|r(z>z.(x)z>i(y)|0 >= ^(A^,(x - + A,^,(x - y)), (14) 
Af (x -y) = %< 0|T(i>,.(x)z>i(l/)|0 >== ^(A^,(x -y)- A^,(x - y)), (15) 



where 
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' {271^ mf -k^ -le 



is the ordinary propagator for a scalar particle with mass rrii. Using the defi- 
nition of 777-1^2 as in Eq. ([l2|) one finds 



' ' (27r)4(mf-A;2-ze)(mi-P-ze) ' ^ ^ 



A^ix) = -ml, / T^T^.^ 7-. ^TT^T^ 7-. r^. (1^ 



d^k e"*''^ 
] {2t^Y {ml - k^ - ie){ml - F - ie) ' 

It is seen that in absence of the (B-L)-violating sneutrino "Majorana"-like 
mass term mj^ = the (B-L)-violating propagator A^^ vanishes while the 
(B-L)-conserving one A? becomes the ordinary propagator of a scalar particle 
with mass rhi = m2 = rfio- 

Majorana neutrino fields can propagate in a virtual state conserving the 
(B-L) quantum number as well as violating it. In the Majorana representation 
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of the neutrino field 



ly^PLiyL + PRM" (19) 

where v'^ — u, the corresponding (B-L)-conserving (S^) and (B-L)-violating 
(S^) propagators can be written as 

S'^ix-y) = i < 0\T{uL{x)uL{y)\0 >^ iPl < 0\T{u{x)u{y)\0 > Pr ={20) 

-tnk^ ^-ik(x-y) 



I 



S^'ix-y) = i<0\T(uL(x)l^{y)\0>^PL<0\T(u(x)P(y)\0>PL^(21) 

_ -ik(x-y) 



d^k 1 



(27r)4m2 - ie 



where = for simphcity. Therefore, the effect of (B-L)-violation orig- 
inating from the neutrino propagator is proportional to the Majorana neu- 
trino mass while a (B-L)-conserving contribution to Oi'f3f3 decay via neutrino 
propagation in the Dirac mode does essentially not depend on the neutrino 
mass and leads to a contribution proportional to the mean neutrino momen- 
tum in a nucleus. The latter is typically of the order of the Fermi momentum 
~ ~ 80MeV. As a result such a contribution, if it exists, is greatly enhanced 
compared to the Majorana mass contribution. 



3 MSSM contribution to Oi/fSP-decay. General 
properties and Effective Lagrangian 

In this section we are considering the general properties of the i?p-conserving 
MSSM contribution to 0z//3/3-decay and derive the corresponding effective La- 
grangian in terms of color-singlet quark charged currents. 

Within our supersymmetric framework there are two sources of lepton 
number violation, the Majorana neutrino mass, m'^, and the (square of the) 
"Majorana" -like sneutrino mass, ml^- Both of them violate - by construction 
- lepton number by two units. Ou^P decay also violates L by two units. The 
OuPP amplitude Roupp will therefore be proportional 

i?0./3/3~m5^OW+ 7^40(2), (22) 

with (9*^*^ representing some matrix elements. Contributions proportional to 
m'^ can be classified by O^^") = 0^^gJ^^ + O^^JjjgY)-> where the first part stands 
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symbolically for the usual mass mechanism diagram, see Fig. 1(a), while the 
second part summarizes all kinds of diagrams involving virtual SUSY particles. 
An example of this type of contribution was given in the introduction, see Fig. 
1(b). Clearly, all diagrams of OuPP decay involving SUSY particles must have 
at least 6 basic vertices. They are thus of higher order compared to Fig. 1(a) 
and can be safely neglected: O^-*^) = O^^Jj^^ + O^^J^gy) — ^\sm)- 

Let us consider i?Oi//3/3 in more details. At the quark level neutrinoless 
double beta decay is induced by the transition of two d quarks into two u 
quarks and two electrons. This process is schematically represented by the 
diagram in Fig. 3(a) encoding all possible contributions to Oufjp-decay. In 
the Rp conserving supersymmetric model it is useful to decompose the basic 
diagram Fig. 3(a) as shown in Fig. 3(b-g). 

Given that the low-energy theory contains a light neutrino there are always 
contributions involving a long-distance interaction component associated with 
the light neutrino exchange. Let us use this fact for a decomposition of the 
quark-lepton Oz^/3/3-cffective vertex in Fig. 3(a) into different parts explicitly 
showing the presence of the neutrino line. Fig. 3(b-d). It is also instructive 
to isolate the SUSY contributions in the form of effective vertices induced by 
heavy SUSY particles exchange. This decomposition is indicated by the small 
black spots corresponding to the short-distance, approximately point-like, ef- 
fective SUSY induced interactions. Fig. 3(c-g). The first diagram. Fig. 3(b), 
corresponds to the conventional standard model Majorana neutrino exchange 
contribution mentioned above. The crossed neutrino line indicates the lepton 
number violating Majorana neutrino propagator, S^. Diagrams Fig. 3(c,d) 
correspond to the neutrino accompanied SUSY contributions. Here, neutrino 
lines are uncrossed and correspond to lepton number conserving neutrino prop- 
agators, S^. (As noted above we neglect SUSY diagrams with lepton number 
violating propagators proportional to the small neutrino mass m^^.) The last 
three diagrams represent the purely supersymmetric contribution and are of 
short-ranged nature. 

This decomposition allows one to apply, if necessary, a Fierz rearrangement 
to the approximately point-like black-spot SUSY vertices in Fig. 3 and rep- 
resent them in the form of a product of color-singlet quark charged currents 
and a leptonic part. Such a representation is crucial for the derivation of the 
0z//9/9-transition operators in the non-relativistic impulse approximation and 
for the subsequent nuclear structure calculations discussed in the sect. 4. 

Assuming the Fierz rearrangement applied to the point-like SUSY vertices 
one can write down the general form of the effective Lagrangian reproducing 



10 



the decomposition in Fig. 3 within 4th order of perturbation theory. It can be 
written in the form: 





■ji ■ eTni)vl + 4 





(23) 



^SUSY 



^SUSY 



+ 



+ 




itW- ■ er,(4)e' 



^SUSY 



'SUSY 



+ 




'jijj ' ^^ij{5)^ ) 



where the SM term is also introduced (see Eq. (pSf ) in Appendix A). 

Color-singlet local diquark operators are defined as 



with a being a color index. The objects rj(fc) and Oi are constructed of Dirac 
gamma matrices as well as derivatives. 

Effective couplings A'^'^^ are dimensionless constants. Different terms are 
scaled out by the characteristic SUSY breaking mass scale msusY with an 
appropriate degree Ui to accommodate correct physical dimension of the cor- 
responding term. As seen from the leptonic part of the effective Lagrangian 
(|23D, the first term conserves lepton number (AL = 0) while the remaining 
five terms violate it by two units (AL = 2). 

In the effective Lagrangian ( pBf ) we neglected possible L-conserving terms 
with the lepton operator structure ~ eFjZ/^,. Their contributions to the Oz//5/3- 
amplitude are strongly suppressed compared to the contributions of the similar 
L- violating terms ~ eTiUl. To see this fact let us have a closer look at the 
corresponding leading order diagrams in Fig. 3(c,d). The bottom parts of 
these diagrams are the SM charged current (SMCC) interactions of the form 
{ifj^PLd){e'j'^PLi^n) ■ Uen, whilc the top parts correspond to the effective SUSY 
vertices. If they are given by the 2nd and 3rd terms of the effective Lagrangian 
Eq. (p3D , the resulting leptonic tensor Csusy can then be written schematically 
as 



Cal=2 ~ nf.PL < 0|T(z/fcZ7„)|0 > PrFc^ • UekU:^ ~ nf^luPRTe' ■ q'/q', (25) 



where q is the neutrino momentum. In the right hand side neutrino masses 
are neglected since m^] << (g), where (g) ~ ^ lOOMeV is the average 
momentum of a neutrino propagating in a nucleus {pp is the nucleon Fermi 



ji = u°'Oid, 



a 



(24) 
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momentum). The mixing matrix elements disappear on the right hand side 
due to the unitarity relation UekU*nSkn = 1- This should be compared with the 
contribution of the possible L-conserving terms ~ eViVi which we neglected 
in the effective Lagrangian in (^). The leptonic tensor in this case takes the 
form 

C-AL=^ ~ eViPL < 0|T(z/fcZ7„)|0 > PLlpC' ■ U^kUen ~ eTapPRe' ■ {m,)/q\ (26) 

where (m^,) = Tn'^M^I^- The same structure appears in the standard neutrino 
mass mechanism with the SMCC at both ends of the virtual neutrino line. 

Comparing Eqs. (^) and (^) one can see that the SUSY contribution 
corresponding to the AL = 2 operators receives from the leptonic sector a 
huge enhancement compared to the contribution of the AL = operators. In 
fact 

C-^L=2lC^L=o ~ PFl{m,) ~ 10' ■ (leV/(m,)) (27) 

For this reason we neglected the SUSY induced AL = operators in the 
effective Lagrangian (PB|). 

Now let us turn from the general consideration to the concrete case of the 
SUSY contribution to Oz//3/3-decay within the MSSM. The following Lagrangian 
terms are relevant to Oz//3/3-quark transitions in Fig. 3(a) 

The Lagrangian terms in the r.h.s. are explicitly given in Appendix A. 

Starting from this Lagrangian one can find 14 dominant diagrams propor- 
tional to which contribute to the Oz//3/3-quark transition in Fig. 3 (a). They 
are listed in Fig. 4. (Note, that in addition to the graphs shown, there exist 
several graphs corresponding simply to an exchange of two of the external mo- 
menta and are not shown for brevity.) As seen, all diagrams in Fig. 4 fall into 5 
classes represented by the last 5 diagrams of the decomposition in Fig.3. The 
supersymmetric part of these diagrams, as discussed above, can be parameter- 
ized by the effective Lagrangian Co^pp given in the general form of Eq. (p3D. 
One can reconstruct a specific form of this Lagrangian in the MSSM comparing 
diagrams in Fig.3 and Fig. 4 and separating the basic SUSY vertices denoted 
in Fig.3 by the black spots representing five different terms of the effective 
Lagrangian C^upp in Eq- (0)- The next step of the derivation is based on 
the standard approximate procedure relying on the fact that all intermediate 
particles involved in evaluation of the effective SUSY vertices are heavy SUSY 
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particles with typical masses of order msusY- As a result these vertices can 
approximately be represented in the form of local operators. The local form of 
the SUSY operators allows one to apply the Fierz rearrangement and to collect 
the quark fields in the color-singlet quark charged currents. Straightforward 
realization of this strategy leads to an effective Lagrangian with the following 
leading order operators violating the lepton number L by 2 units 

n^ /c^ W~W~ 

Cal=2 = -{v^X + V^ww + V%^w)— e(l + 75)e^+ (29) 

^SUSY 

''''SUSY 

The color-singlet quark charged currents are defined as usual 

j^, = cos^, ^7^(1-75)^^. (30) 

Note that since we take only the leading order contributions into account in 
Eq. (^91) not all possible terms of the decomposition Eq. (^) are retained in 
Eq. (|29|). Naively one might have expected that the diagrams in Fig. 4 (and 
the corresponding terms in Eq. (|23|) ) are ordered with respect to decreasing 
importance, since a larger number of heavy sparticles in the loops results in 
larger (loop) suppression factors. However, the exphcit calculation shows that 
this is not the case. Terms corresponding in structure to the 2nd, 3rd and 
the 5th terms in Eq. (p3D are suppressed by the helicity structure of the basic 
MSSM interactions and/or typical suppression factors of order rrie/msusY from 
left-right sfermion mixing or higgsino-like interactions. 

Retaining only leading contributions to the operator structures in Eq. (p9|), 
the dimensionless lepton number violating parameters take the form 



72 "^msvsY' il ^msusY>^msvsY> 



( "^^ )g{mg,m ±,m ±, 
^msusY^ 



Vgu 



^(^)' E ViiVji{^^)r(m„ m^.,m^.), (32) 

72 ^msusY^ Xj ^^susY^ 

4 ^rnsusY^ 
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TTl lb 777, i 

'^msusY^^msusy'^msusY^ 
= -r( —) Y.^(^e,m^^,m ±)eLSe)VjiX (34) 

4 ^TflsuSY ' 



o«(^) + o,S(^) 

4 ^rUsuSY^ i ^^SUSY^ 

The dimensionless loop factors T {nii) , Q {mi) , J {mi) and X{mi) are given in 
Appendix B. They depend on the sparticle masses mj in the corresponding 
loop. Recall that g and Qs are the SU{2)l and S'f/(3)c coupling constants. 
Further definitions on couplings and mixing parameters can be found in Ap- 
pendix A. 

The next step of the calculation deals with reformulating the problem in 
terms of nucleon degrees of freedom instead of quark ones. This is relevant for 
the nuclear structure part of the calculations. 



4 From quark to nuclear level 

So far the discussion has focussed on particle physics aspects, deriving the low- 
energy effective Lagrangian in Eq. ( P^D formulated in terms of quark fields. 
However our goal is the calculation of the amplitude Ro^pp of Oz//5/5-decay which 
is a nuclear process proceeding not at the level of quark degrees of freedom 
but at the level of nucleon ones. Formally one can write down 

7^o./3/3 = <{A,Z + 2),2e-\S-l\{A,Z)>= (36) 
= < {A,Z + 2),2e-\Texp[i J d''xCo,pf3{x)]\{A, Z) > 

where the effective Lagrangian £0^/3/3 = l^wff + '^al=2 is given by Eqs. 
(pSD, (^). The nuclear structure is involved via the initial (A,Z) and the final 
(A, Z-l-2) nuclear states having the same atomic weight A, but different electric 
charges Z and Z-l-2. The standard framework for the calculation of this nuclear 
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matrix element is the non-relativistic impulse approximation(NRIA). It implies 
the substitution of the quark current j^Av in the effective Lagrangian £oi^/3/3 
in Eq. (PB|) by the non-relativistic nucleon current 



NRIA jf, 

J AV * ^ • 

The latter is an incoherent sum over individual nucleon currents of a nucleus 
and is given by the formula |15 



J'^(x) 



Here fy ~ 1, /a ~ 1.261, g^^ is the metric tensor, is the isospin raising 
operator, is the position of the ith nucleon, the superscript k stands for the 
spatial component. Ci and Dj are the well-known scalar and the vector nuclear 
recoil terms given by B] 



D, 



2m p 
1 

2m p 



IP* + Pi 



fA 



[Pi + Pi) + ^(1 - 2mp^)qi X or. 

Jv 



(38) 
(39) 



Here {pi,Ei) and (p-,-E'-) are initial and final 3-momentum and energy of 
the ith nucleon and the 3-momentum transfer is = Pi — p^. The nucleon 
couplings obey the relations 



fw/fv = -(/ip - lJ^n)/i'2mp) ^ 3.7/(2mp), /p//a = 2mp/m. 



2 



(40) 



where mT^ is the pion mass and /ip(n) is the proton (neutron) magnetic moment. 

Since we are interested in the dominant contributions only, recoil terms 
will be neglected in the rest of this paper and have been given above for 
completeness. 

The exponential factor in Eq. (|37| ) is introduced instead of the local delta 
function in order to take into account the finite nucleon size. It is the Fourier 
transform of the nucleon form factor F(q ^) in the conventional parameteriza- 
tion a dipole form 



F(q^)= 1 + ^ 



m. 



(41) 
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with tua = 0.85GeV. The finite nucleoli size effects are known to be 



important for the short-distance contributions to the Oi//3/3-amphtude such as 
those corresponding to the dominant terms in the effective Lagrangian 

Now, starting from Eq. (p6D, it is straightforward to calculate the Ou^jS- 
amplitude Rquis/b within the non-relativistic impulse approximation. The final 
result for the 0"'' — > 0"*" transition amplitude can be written as follows 

^SUSY 

RouppiO^ ^ 0+) = V2 ■ e(l + 75)e^ < F\n'^'^\I > . (42) 

"^SUSY 



The normalization factor is 

47r Ro 
nipnie /a 



Cou = z— ^- (43) 



Here, Rq is the nuclear radius, mp and rrie are the proton and the electron 
masses. 

The effective lepton number violating parameter is defined as 

V'"'" = + + + (44) 

In Eq. we have introduced the transition operator Q^^^"^ in order to 
separate the nuclear physics part of the calculation from the particle physics 
one. Having the transition operator one can calculate the corresponding nu- 
clear matrix element for any Oi//5/5-decaying candidate isotope within any spe- 
cific model of nuclear structure. From now on we define 



This nuclear matrix element is found to be equal to those for heavy neutrino 
exchange [^] 

j^susY ^ |_^^^^ _ Mgt,n}, (46) 

where 

2 

Mgt,n = (^^) {F\nGT,NU) (47) 

where 

{F\naT,N\I) = {F\T.^?r^^^^^ " f") M^aW (48) 
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and 



where 



M 



'fv^2, 
fA 



{F\n,Mi) = {f\y: 



Tl'T 



+ 




FNixAm. 



Here, F/v(xa) is the short-ranged potential 



Fn{xa) 



1 



,2M 



(49) 



(50) 



(51) 



(27r)3 {m\ + 

mATij. This potential takes into account the finite nucleon size, see 



with xa 

Eq. (41), its analytic solution is given by the formula 

X 



Fn{x) 



48 



(3 + 3a; + x^)e" 



(52) 



The above definitions are general in the sense that one can apply nuclear 
wave functions of any nuclear structure model for their calculation. For the 
following analysis, on the other hand, numerical values for the matrix elements 
are needed. 

In our numerical analysis we will use the following value for the Oz//3/3 decay 

of ^^Ge rri 



M 



SUSY 



289. 



(53) 



This value is based on a pn-QRPA model [1^, which has been discussed al- 
ready several times in the literature |jl8[ and has been applied previously to 
calculations of the R-parity violating contributions to Ou^P decay H, |lOl, as 
well as to Oz//3/3 decay in left-right symmetric models |T^. We will therefore 
not repeat the details of the calculation here, and refer for brevity to [0 . Un- 
certainties associated with the pn-QRPA have been discussed in [Q for J^f,n 
and ^AGT n- 



5 Ojy/3/3 constraints on (B-L)-violating Sneutrino 
Mass 



The OuPP-decaj amplitude given in Eq. (|42|) leads to the following half-life 
formula 

2 



[r°;f (0+ - 0+)] 



4ml 
'-^01 7^ 



7] 



SUSY 



-M 



SUSY 



m 



SUSY 



(54) 
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Here Gqi is the standard phase space factor tabulated for various nuclei in . 

Eq. (^) takes into account only the contributions from sneutrino exchange. 
There might be other contributions which we assume not to cancel the SUSY 
contribution. If there is no unnatural fine tuning between different contribu- 
tions we may retain only the SUSY one in deriving upper bounds for the lepton 
number violating parameters. 

The most stringent experimental lower limit on 0i//3/9-decay has been ob- 
tained for ^^Ge 

T°f ^-""P(0+ ^ 0+) > 1.0 X lO^Vars 90% c.l. (55) 

Combining this bound with Eq. (|5^) and the numerical value of the nuclear 
matrix element Ai^^^"^ given in Eq. we get the following constraint on 

the effective MSSM parameter 

Since we are interested in deriving constraints on the (B-L)-violating sneutrino 
mass rfiM from fj^^^^ , we will adopt the following simplifying assumptions. 
Assume all SUSY particle masses to be equal to the effective SUSY breaking 
scale msusY introduced in Eq. (|23|) and consider two limiting cases for the 
lightest neutralino x composition. In the first case it is assumed to be a pure 
B-ino as suggested by the SUSY solution of the dark matter problem |]19|, in 
the second a pure Higgsino. These two cases can be understood as extreme 
cases, and actual values for ttlm for other choices of the neutralino composition 
should therefore lie in between the two extreme values given below. In the 
Higgsino case essentially only the last three graphs in Fig.4 with gluino lines 
contribute to 0z//3/3-decay. As seen from Eq. (^) the corresponding lepton 
number violating parameter T]q does not depend on the neutralino composition 
and survives in this limiting case. With the currently accepted values of 
the gauge coupling constants and W-boson mass we derive 

6 Neutrino mass constraints 

As already mentioned, the sneutrino contributes to the Majorana neutrino 
mass m'^j at the 1-loop level via the (B-L)-violating propagator Eq. 
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proportional to the sneutrino (B-L)-violating mass parameter The corre- 
sponding diagram given in Fig. 2(f) gives rise to an induced Majorana neutrino 
mass 6m'^. Thus, in the presence of a non-zero tree-level contribution m^J^^^'^^ 
the total neutrino mass is 

mlj = Tn^'''^'^ + 6nf. (59) 

As seen from Eq. (|T^ the (B-L)-violating sneutrino propagator in momentum 
space is a rapidly decreasing function of momentum q. In the ultraviolet limit 
it behaves as A^^{q) ~ unlike an ordinary scalar propagator decreasing 
only as ~ As a result, the 1-loop diagram in Fig. 2(f) leads to a finite 

loop integral. Hence, the sneutrino induced Majorana neutrino mass dnf is a 
calculable object. It is given by the formula 

4 

= ^ ef^m^,m^^^j/(m^, m^, m^J (60) 

k=l 

where the subscript i stands for generation. The neutralino-neutrino-sneutrino 
coupling is e^^ = tan 6*^/ A/la — ■^k2- The loop integral is defined as 

j,~2 ~2 N_ ■ { d^<l 1 1 1 

J [nil ~ Q ) ("^2 ~ Q ) i^xk ~ ^ / 

For an approximate numerical estimation we take all superpartner masses equal 
to the common mass scale of supersymmetry breaking m^^. ~ rhi ~ m2 ~ 
msusY- Then in this approximation one gets for the lightest neutralino x 
contribution the following constraint on the sneutrino mass splitting parameter 

"^^«^|tan^H.Arn-Ar,,|llOOG^j ^^^^ ^''^ 

Here m^j^^.^ < ml^^ are the best laboratory limits on the neutrino masses which 
can be summarized as lEOl mff = 15eV, mf,^^ = 170KeV, mff^ = 23MeV. 
Eq. (|6^ ) assumes that there is no significant cancellation in Eq. (^91) between 
the tree and 1-loop level contributions. 

If the neutralino is B-ino dominant, than we derive the following limits on 
the sneutrino mass splitting parameter 

rfiMie) < 120MeV, rfiMi^) < 13GeV, mM{r) < 149GeV. (63) 
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Thus, for the second and third generation sneutrinos large sphttings are not 
excluded by experimental data. 

An interesting question to ask is whether there are certain loopholes in 
the constraints on the sneutrino mass splitting derived from the experimen- 
tal upper limits on neutrino masses. It could happen, for instance that the 
lightest neutralino is higgsino dominated, in which case one would expect that 
the bound (|6^) might have to be relaxed considerably. However, one should 
remember that all neutralino states contribute to (0), so that even if there 
is no constraint from the lightest neutralino the other mass eigenstates will 
provide a finite contribution to the neutrino mass. 

In order to investigate this question a little bit more quantitatively we did 
a numerical scan of the SUSY parameter space, calculating upper bounds on 
ifiM taking into account all four neutralino states. In this case instead of Eq. 
we have 

~ < 1.7-10-^GeV (rnsusY_\ f^V^' 

"""'^'^ - {Et=i{t^n9wNa-N,,yy,Cs{x,yi)y/^ UoOGeV^ 1^ leV j " ^ ' 

Here x = rhn /msusY and yi = m^JmsusY- C^{x,yi) takes into account the 
fact that the masses of the particles in the loop are no longer taken to be equal. 
In the limit where friM <^ ttld the 1-loop integral C^{x,yi) is given by 

r i X + ^^(Log(l/^) - Log(x^) - 1) 

Cz[x,y) = 2 — — (65) 

C3{x,y) is normalized such that it approaches 1 in the limit where x and y 
approach 1. We let the parameters of the neutralino mass matrix vary from 
(0 — 1000) GeV for fi and M2, tan/5 from (1 — 50), for both positive and 
negative fi. The unification condition Mi = (5/3) tan^ 6wM2 has been assumed 
in this calculation. We required the lightest mass eigenstate to be heavier 
than 20 GeV, motivated by the LEP measurements. About 10* solutions 
were calculated. From these we calculated the "average constraint" and an 
"absolute" upper bound. These are: 

< 60(125) f^j MeV. (66) 

on average ("absolute"), if m^i ~ msusY = 100 GeV is assumed. These 
numbers are about a factor of 2 (4) less stringent than taking only the lightest 
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neutralino (being bino) fixed at 100 GeV. This simply reflects the fact that 
within the above-mentioned parameter ranges many solutions exist where even 
the lightest neutralino mass state can be considerably heavier than 100 GeV. 
On the other hand, it seems that within the typical range of SUSY parameters, 
the constraint on rfiM is essentially "stable" and has to be taken seriously. 
(This conclusion remains unchanged even if we drop the uniflcation assumption 
on Ml, although the bounds might have to be slightly relaxed in some cases.) 

Note that, in principle, more stringent limits on rfiM than in Eq. (|63D 
could be derived from the upper bounds on the neutrino mass given by non- 
observation of Oi//5/3 decay. However, in this case the situation is more complex, 
since 0^(3(3 decay measures an effective neutrino mass {m'^) = J2' Uljirij, where 
Uej are mixing coefficients connecting the weak and the mass eigenstate basis 
for neutrinos. Thus limits on rfiM derived from the neutrino mass limit of Qv(3(3 
decay will also depend on assumptions on neutrino mixing. If one assumes for 
simplicity Uej ~ ^ei one could derive mM(e) < 22 MeV from the data on ''^Ge 

i- 



7 Conclusion 

In summary, we have proven a low-energy theorem for weak scale softly broken 
supersymmetry relating the (B-L)-violating mass terms of the neutrino and 
the sneutrino as well as the amplitude of neutrinoless double beta decay. This 
theorem can be regarded as a supersymmetric generalization of the well-known 
theorem |]14| relating only the neutrino Majorana mass and the neutrinoless 



double beta decay amplitude. 

According to Eq. ([T2|) the parameter rfiM describes a splitting in the 
sneutrino mass spectrum. This splitting leads to mixing in the sneutrino- 
antisneutrino {v — v^) system and to the effect of lepton number violating 
p-p^ oscillations f^, [^]. 

The mass splitting parameter friM is constrained by the experimental data 
on neutrinoless double beta decay Oi//3/3 and the neutrino mass discussed in the 
present paper. The neutrino mass constraint on rfiM is found to be more strin- 
gent then the direct Oi//5/3-decay constraint. This is opposite to the conclusion 
reached for R-parity violating SUSY, for which the direct double beta decay 
constraints have been found to be more stringent than those derivable from 
the neutrino mass However, in contrast to the neutrino mass limits, 

the corresponding constraint on niM from neutrinoless double beta decay is 
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completely independent of assumptions about neutralino masses and mixings. 

The constraints derived here, nevertheless, leave quite some room for ac- 
celerator searches for sneutrino mediated (B-L)-violating effects for the second 
and third generation. The sneutrino mass splitting parameter rfiM might be 



searched for at future colliders such as the NLC or a first muon collider 21 



2^. Probably, dedicated searches for Majorana sneutrinos have a chance of 
detecting positive signal, within the above discussed low-energy limits. For 
z/g, on the other hand, these limits seem to be too stringent and accelerator 
experiments should not be expected to give positive signals. 
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8 Appendix A. Supersymmetric Lagrangian terms 
contributing to the Oupp-decay amplitude 

In the presence of the (B-L)-violating (s)neutrino masses, given in Eq. (0), 
Oi//3/3-decay is triggered by the following terms of the MSSM Lagrangian 

^MSSM = ^Wff + ^Wff + ^X+ff ~^ ^Xff ~^ ~^ ^Wx+X ^^'^) 

The individual terms can be found in the standard sources like ref. [0, p3| . 
Let us list them explicitly. 

a) Gauge boson-fermion-fermion 

This is nothing but the usual standard model charged-current Lagrangian, 



r - - - ^ 

^Wff - ^ 



W;l{uLj>'dL + plI^cl) + h.c.J (68) 

= ^^Wqq + C-WTl- (69) 

b) Gauge boson-sfermion-sfermion 
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Only the charged- current part of this type of the MSSM interactions is of 
interest in OujSP decay: 

^Wff — ^Wqq + ^wli i^^) 
= -^W^iul'i'd,) - ^W^ii>0'e,) + h.c. (71) 

Note that d'^ is defined by = — . 
c) Chargino-fermion-sfermion 



C^^ff- = Cl^-ULxtdL + Ci^-dLXr^L + C^LR-^LxtdR+ (72) 

+ C^L ■ URxtdi + Cj^L ■ dRX'UL + Ci^^ ■ dLXiUR + 
+ ClL-^LxteL + ClL-eLXt^L + ^-c, 

where the following shorthands have been defined: 
Cll = -gUii, Cil = -gVa, 

_ gmd ^d _ gmu 

y/ 2mwcosp y imwsmp 

Coefficients C^l, Crr are fermion-sfermion couplings to the gaugino compo- 
nent of the chargino while Clr,Crl describe couplings to the Higgsino com- 
ponent. The latter are proportional to the fermion mass and, therefore, can 
be neglected in the fermion-sfermion sector as is done in the present paper. 

The chargino mixing matrices Uij and Vij are defined from the diagonal- 
ization of the chargino mass matrix M^± 

U* ■ M^±V^ = Diag{m^±}, (73) 

For details see Ref. |23|1 . 

d) Neutralino- fermion-sfermion 
The neutralino interaction can be written as 

^xff = [4(i)^LXi^L + 4i{i)^RXi'4^R + (74) 

+ 4R(i)i'LXii'R + 4iL{i)^RXii^L\ - h.c. 
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where iIjl = UL,dL,eL,i^L and their scalar superpartners ipL = UL,dL,eL,i^L- 
The corresponding chiral coefficients are 



4a) = Qm^riOwMi, (75) 



(76) 



mwsin(3 

^^Ar,4, (77) 
AT-, (78) 



4«., = ^^-V,3. (79) 



my/cosp 

Coefficients Nij are elements of the orthogonal neutralino mixing matrix 
diagonalizing the neutralino mass matrix M^. In the MSSM the neu- 
tralino mass matrix is identical to the MSSM one [|l^] and in the basis of 
fields (S, W^, Hi H^) has the form: 



0) 



/Ml -Mzs^cp Mzs^sp \ 

M2 Mzc^c^ -Mzcy^s^ 

-MzS^rC/3 MzC^Cp -jl 

V Mzs^sp -Mzc^sp -fi y 

where Cy^ = cos 9w, s^r = sinOw, ty^r = tan6'v^', sp = sin/5, = cos/?. The 
angle P is defined as tan (3 =< H^ > / < H^ >. Here < H^ > and < H^ > 
are vacuum expectation values of the neutral components H2 and H^ of the 
Higgs doublet fields with weak hypercharges Y[H2) = +1 and Y{H^) = — 1, 
respectively. The mass parameters Mi,M2 are the soft SU{2)l and U{1)y 
gaugino masses. In grand unification scenarios they are related to each other 
as follows 

5 

Mi = -ian^Ow ■ M2 (81) 
3 

By diagonalizing the mass matrix (^) one can obtain four neutralinos Xi with 
masses m^. and the field content 

X^ = UilB + M^2W^ + Mgijf + Ui^Hl (82) 

Recall again that we use notations W^, B for neutral SU{2)l x U{1) gauginos 
and Hi H2 for higgsinos which are the superpartners of the two neutral Higgs 
boson fields H^ and H^ 
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We apply a diagonalization by means of a real orthogonal matrix M. There- 
fore the coefficients Mij are real and masses m^- are either positive or negative. 
The sign of the mass coincides with the CP-parity of the corresponding neu- 
tralino mass eigenstate Xi- If necessary, a negative mass can be always made 



positive by a redefinition of the neutralino field Xi- It leads to a redefinition 
of the relevant mixing coefficients Mij — >■ i ■ Mij. 

e) Gluino-s quark- quark 

The gluino interaction is given by 

= -V2g,^ - rR-g^^^fR) + h.c, (83) 

Here A*^"-* are 3x3 Gell-Mann matrices (a = 1, ...,8). Superscripts a,P in eq. 
(|83D are color indices. 



f) Gauge boson- chargino-neutralino 
In the notation of ref. ||T3| 



where 



= 9W-x^Y{O^^PL + OfjPR)xt + h.c, (84) 



Ot, = -^MuV;,+M..2V;, (85) 



1 

Further details and conventions on the definitions used can be found in the 



OS = +^m:,u,,+m:,u,,. (86) 



paper by Haber and Kane |13|. 



9 Appendix B. Box Integrals 

In this appendix some relevant formulas for the calculation of the loop integrals 
are summarized. There are four kinds of integrals encountered in the Eqs. 
5TD-(^: 



g{mi,m2,m3) = -i j -^^-^ (87) 



X 



10 

SUSY 



{nil ~ k'^){m2 — k'^){rn?-r — k'^Yimi — k'^){m2 — k'^){m1 — k'^) 
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J^(mi, 1112,1713) = y (88) 
(mf — k^){m2 — A;^)(m| — k'^y{m\ — k^){m2 — k^){ml — k^) 



X(mi,m2,m3) = ~^ / (89) 



(27r) 



X 



''^SUSY 



(mf — A;^)(m2 — A;^)(mf — A;^)(m2 — A;^)(m3 — A;^) 



J{mi,m2,m3) ^ -'^ J ^^^-^ 

^ ''^SUSY^'^ 

{ml — A;^)(m| — A;^)(mf — A;^)(m2 — A;^)(m| — k'^) 

All four integrals are finite and can be calculated analytically using stan- 
dard methods. Simple solutions can be found for the case when the masses of 
all particles in the loops are assumed to be equal to mgusY- In this limit one 
finds 

G{msusY) = (4807r2)-\ r{msusY) ^ {^^QT^''y\ (91) 
AmsusY) = J(m5c/5y) = (1927r')"'. 
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Figure Captions 

Fig.l.: (a) the conventional Majorana neutrino mass mechanism of OuPP- 
decay; (b) an example of a Rp conserving SUSY contribution to the OuPP- 
decay. 

Fig. 2.: Lowest order perturbation theory diagrams representing the relation 
between the neutrino Majorana mass m'^, the " Majorana" -hke (B-L)-violating 
sncutrino mass ffiM and the amplitude of neutrinoless double beta decay Rovpf}- 
(a) the neutrino and (b) an example of sneutrino contribution to the Oi'f3f3- 
decay amplitude Rovpp- 0i//3/3- vertex contribution to (c) the neutrino Majorana 
mass and (d) to the "Majorana" -like sneutrino mass; (e) neutrino contribution 
to the sneutrino "Majorana" -hke mass and (f) sneutrino contribution to the 
neutrino Majorana mass. Crossed (s)neutrino lines correspond to the (B-L)- 
violating propagators. 

Fig.3.: The decomposition of the Oi//3/3-effective vertex, (a) to the left: 
The Oi//3/3-effective vertex. To the right the decomposition: (b-d) first line: 
neutrino- accompanied contributions to {)vj3j3 decay, (e-g) second line: purely 
supersymmetric contributions without neutrinos. 

Fig. 4.: Rp conserving MSSM contributions to the Oi//3/3-decay amplitude. 
Leading order diagrams, see also the text. 
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Figure l.(a) 
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